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ABSTRACT 

The steady-state structure of a disc with a corona is analyzed when the vertical com- 
ponent of the gravitational force due to the self-gravity of the disc is considered. For 
the energy exchange between the disc and the corona, we assume a fraction / of the 
dissipated energy inside the accretion disc is transported to the corona via the mag- 
netic tubes. Analytical solutions corresponding to a prescription for / (in which this 
parameter directly depends on the ratio of the gas pressure to the total pressure) or 
free / are presented and their physical properties are studied in detail. We show that 
the existence of the corona not only decreases the temperature of the disc, but also 
increases the surface density.The vertical component of the gravitational force due to 
the self-gravity of the disc decreases the self-gravitating radius and the mass of the 
fragments at this radius. However, as more energy is transported from the disc to 
the corona, the effect of the vertical component of the gravitational force due to the 
self-gravity of the disc on the self-gravitating radius becomes weaker, though the mass 
of the fragments is reduced irrespective of the amount of the energy exchange from 
the disc to the corona. 
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1 INTRODUCTION 

Accretion plays a crucial role in the processes of energy lib- 
eration and mass accumulation near the cores of compact as- 
tronomical objects (Bertout 1989; Robinson 1976; Hunter, 
Ball & Gottesmanet 1984; Kohler 1995), and in particu- 
lar for Active Galactic Nuclei (AGN) (Shields 1978 ;Os- 
terbrock 1993;Szuszkiewicz 2001; Czerny 2007). Self-gravity 
has global consequences on the disc shape, influencing the 
location of its inner and outer edges, as well as the disc geo- 
metrical thickness (Karas, Hure & Semerak 2004). The role 
of self-gravity has been separated into that from its vertical 
component (e.g., Fukue & Sakamoto 1992; Hure 1998) or 
from its radial component (Lu, Yang & Wu 1997), and from 
combination of two (Yang & Liu 1990; Bertin & Lodato 
1999). The vertical structure of accretion discs of galactic 
nuclei is affected by their own self-gravity and the gravita- 
tional force of the central object, especially in the middle 
and outer regions, where the vertical component of the disc 
gravitational field exceeds the corresponding vertical com- 
ponent due to the central object. 
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In addition to the role of self-gravity, one of the main 
open issues in the physics of black hole accretion discs is the 
relationship between the disc MRI-driven turbulent viscos- 
ity and the generation of a hot corona, which is postulated 
in order to explain the observed X-ray emission (Liang & 
Price 1977; Galeev, Rosner & Vaiana 1979; Blackman & 
Field 2000; Kuncic & Bicknell 2004). Furthermore, theoret- 
ical simulations successfully reproduce the observed spec- 
tra by assuming a certain fraction of gravitational energy is 
released as hot gas (e.g., Haardt& Maraschi 1991; Svens- 
son & Zdziarski 1994; Kawaguchi, Shimura & Minishige 
2001) and/or a certain spatial distribution of the hot gas 
(e.g., Merloni & Fabian 2001b). Svensson & Zdziarski (1994) 
studied a disc-corona system, in which a major fraction of 
the power released when the cold matter accretes is trans- 
ported to and dissipated in the corona. They showed that the 
steady-state structure of the disc significantly modified due 
to the existence of a corona. Miller & Stone (2000) showed 
that the primary saturation mechanism of the MRI is local 
dissipation, about twenty five percent of magnetic energy 
generated by MRI within two scale heights escapes because 
of buoyancy, producing a strongly magnetized corona above 
the disc. It was shown by Merloni & Fabian (2001a) that the 
inferred thermal energy content of the corona, in all black 
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hole systems, is far too low to explain their observed hard 
X-ray luminosities, unless either the size of the corona is at 
least of the order of 10 3 Schwarzschild radii, or the corona 
itself is in fact a reservoir, where the energy is mainly stored 
in the form of a magnetic field generated by a sheared rota- 
tor. 

Following these considerations Khajenabi & Shadmehri 
(2007) (hereafter KS) studied the gravitational instability 
of an accretion disc with a corona. They showed that the 
disc becomes more gravitationally unstable due to the ex- 
istence of the corona. Because of the corona, not only does 
the self-gravitating radius decrease, but also the mass of the 
fragments at this radius increases (KS). However, KS did 
not take into account the vertical component of the gravi- 
tational force due to the self-gravity of the disc. Also, their 
analysis was based just on the Toomre parameter criteria. 
Actually, theory and simulations (Gammie 2001; Rice et al. 
2003, 2005; Lodato & Rice 2004; Mejia et al. 2005) show 
that two conditions should be fulfilled for accretion discs to 
fragment gravitationally. First, the accretion disc must be 
massive enough so that gravity can overcome thermal pres- 
sure and centrifugal support, which implies a value of less 
than unity for the Toomre parameter (Toomre 1964). Sec- 
ond, the disc must cool fast enough for the compressional 
energy provided by the collapse to be radiated away. So, ac- 
cording to the numerical simulations (Gammie 2001; Rice 
et al. 2003, 2005; Lodato & Rice 2004; Mejia et al. 2005), 
the cooling time-scale must be of the order of the dynamical 
time-scale. For an accretion disc with a corona, we believed 
that the corona provides another mechanism for the cool- 
ing of the disc and so, existence of the corona decreases the 
cooling time-scale. 

In this paper, we extend the analysis of KS in different 
directions: (i) Toomre parameter and the cooling time-scale 
conditions are considered as criteria for the fragmentation 
of the disc; (ii) the vertical component of the gravitational 
force due to the self-gravity of the disc is considered; (iii) for 
the fraction of the dissipated energy inside the disc that is 
transported to the corona, we consider two cases, in the first 
case, this fraction depends on the physical variables of the 
system (Merloni & Nayakshin 2006); and in the other case, 
this fraction is assumed to be constant and we consider it as 
a free parameter of the model. In the next section, the basic 
assumptions and the equations are presented. We analyse 
properties of the solutions in Section 3. The gravitational 
stability of our solutions is studied in detail in Section 4 
. By studying the gravitational stability of such a system, 
we address the effect of self-gravity on disc stability and 
fragmentation. We conclude with a summary of the results 
and a discussion in the final section. 



effects are only important at small radii. Thus, the rotational 
angular velocity o f the disc is approximately Keplerian, i.e. 
n K = \/ GM/R 3 . Although we neglect the radial compo- 
nent of the self-gravity of the disc,the self-gravity of the disc 
in the vertical direction that modifies the hydrostatic equi- 
librium of the disc is included in our model. We take it into 
account with the extra assumption that the accretion disc 
is homogeneous and extends to infinity, as done in several 
other models (Paczynski 1978; Kozlowski, Wiita & Paczyn- 
ski 1979; Sakimoto & Coroniti 1981; Shore & White 1982; 
Cannizzo & Reiff 1992; Hure et al. 1994a). The weight of 
the corona is negligible as has been shown by Svensson & 
Zdziarski (1994). So, the vertical hydrostatic equilibrium of 
the disc implies 
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which measures, along the axis of rotation, the accelera- 
tion due to self-gravity as a fraction of the central gravi- 
tational acceleration. With this definition, equation (1) can 
be rewritten as 
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In what follows, when we discuss solutions without self- 
gravity of the disc, we mean solutions that are obtained 
by neglecting £ in the above equation (see KS). Also, by 
solutions with self-gravity of the disc, we mean solutions 
that are obtained by considering ( in the above equation. We 
obtain solutions with self-gravity in this paper and illustrate 
the possible effect of the self-gravity by making comparisons 
to solutions without self-gravity. 

The basic remaining equations of our model are similar 
to KS. However, we briefly mention the equations and basic 
assumptions for completeness. We consider a more general 
prescription for the viscous stresses r r ^ (Taam & Lin 1984; 
Watarai & Mineshige 2003; Merloni & Nayakshin 2006): 
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where ao and < fi < 2 are constants and p is the sum 
of the gas and radiation pressures. Phenomenological mod- 
els generally assume that at each radius, a fraction / of 
the accretion energy is released in the reconnected magnetic 
corona. Assuming that in MRI-turbulence discs such a frac- 
tion / of the binding energy is transported from large to 
small depths by Poynting flux, Merloni & Nayakshin (2006) 
estimated the fraction / as 



2 GENERAL FORMULATION 

As a result of differential rotation the gas is heated, loses an- 
gular momentum and a slow (subsonic) radial drift results. 
In a steady state the mass accretion rate is constant, and 
the local viscous energy release is balanced by the radiative 
cooling. In the direction perpendicular to the plane of the 
disc, hydrostatic equilibrium is assumed. The rotation curve 
is dominated by a Newtonian point mass M, as relativistic 



(5) 



where (3 is the ratio of gas pressure to the total pressure 
(The MN prescription hereafter). With expressions (4) and 
(5) for the viscous stresses and dissipated power respectively, 
we can construct the basic equations describing the disc. The 
azimuthal component of the equation of motion gives 
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4crf2 K , 



where J(R) = I — ^/R^JR and R lrl denotes the inner bound- 
ary of disc. Since we are interested in the regions of the disc 
with radii much larger than R 1D , we have J « 1, but we 
keep the more general expression in what follows. Energy 
conservation implies that 
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Finally, with the vertical transport of heat dominated by 
radiative diffusion we have a relation between the midplane 
and surface temperatures given by 



T = (^E) 1/4 T cff 



(8) 



where n is the opacity coefficient. 

Equations (3), (6), (7) and (8) enable us to find p and 
T and p as functions of R with f3 and C, the critical input 
parameters. Thus, 
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There are two algebraic equations for /3 and £ as follows 
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where /i m is the mean particle mass in units of the hydro- 
gen atom mass, mn- The other constants have their usual 
meaning. These two equations give us £ as follows 
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Now equation (13) enables us to find T, P and p as functions 
of J?, /3 and the other input parameters in the form 
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In order to study the behavior of our solutions, it is 
more convenient to introduce dimensionless variables. For 
the central mass M, we introduce Ms = M/(1O 8 M ) and for 
the radial distance R, we have r-j, — i?/(10 3 Rs), where Rs = 
2GM/c 2 is the Schwarzschild radius. The mass accretion 
rate can be written as 

1-e A-kGM _ l E L E 



M 



(17) 



where Z E = L/L E is the dimensionless disc luminosity rel- 
ative to the Eddington limit, e = L/(Mc 2 ) is the radiative 
efficiency and K e . s . ~ 0.04 m 2 kg _1 is the electron opacity. 
In our analysis, we will use the nondimensional factor Ze/e 
as a free parameter so that by changing this parameter we 
can consider appropriate values of the accretion rate. How- 
ever, some authors introduce different forms for the accre- 
tion rate. For example, Nayakshin & Cuadra (2005) who 
studied gravitational stability of the Galactic Center, intro- 
duced M = (rh/e)(L E /c 2 ) with e w 0.06 and m = 0.03 to 1. 
These values, which are appropriate for the Galactic Cen- 
ter, correspond to Ze/e ~ 0.5 to 16.6 in our notation. Also, 
for a central mass with mass M — 10 8 Mq, Goodman & 
Tan (2004) proposed Z E /e = 10. Thus, in our analysis, the 
corresponding range of values for Z E /e is from 1 to 10. 
Our solutions for T, P and p become 
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We can also calculate the surface density as 

E = 3.292 x 10 5 a - 4/5 ^ 1/5 M 8 1/5 (^) 3/5 J 3/5 r 3 - 3/5 
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and the ratio /3 is obtained from nondimensional form of 
equations (11) and (13), i.e. 
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C = 0.134 x a ( 7 3/5 « 3/5 M 8 7/5 (^) 6/5 J 6/5 r 3 3/10 
x/3 (i6-3 M )/io (1 _ ^-i^ _ y / 2 a f3^ 2 y i/5 (23) 
0.216 x a; /5 ^ 5 M 8 1/5 (^) 8 / 5 J 8 / 5 r 3 - 21/1() 

x/3 (s+m)/io ( - 1 _ ^-2(1 _ v /2a /3' l/2 ) 9/5 - 1 - C = (24) 

where k = /t//t e .s. and we assumed /i m = 0.6. 

Equations (18), (19), (20) and (21) describe the struc- 
ture of a self-gravitating disc with a dissipative corona. How- 
ever, the physical variables depend not only on the radial 
distance but also on (5 and £ which can be calculated at 
each radius from algebraic equations (23) and (24). In the 
next section we will analyze our solution, in particular the 
gravitational stability of the disc. 



3 ANALYSIS 

We can solve algebraic equations (23) and (24) numerically 
at each radius to obtain the ratios (5 and £. Having these pa- 
rameters, physical variables are obtained based on our ana- 
lytical solutions. Since we are interested in gravitational sta- 
bility of the disc, we can approximate J{R) ~ 1 for _R 3> _Ri n . 
Also we fix the central mass Ms = 1, the opacity k = 1 and 
the viscosity exponent fi = 1, but vary the other input pa- 
rameters. However, for the transport of energy from the disc 
to the corona (i.e. the fraction / in our approach), we may 
either consider the MN prescription (in which / is not con- 
stant) or assume that the fraction / is constant and use 
it as a free parameter of the model. We study properties 
of the solutions in both cases in the next two subsections. 
Interestingly, typical behaviors of the solutions are similar 
qualitatively, though precise values of the physical variables 
at each radius depend on the mathematical prescription of 
the fraction /. The possible effects of the self-gravity of the 
disc in both forms of prescription for the fraction / are il- 
lustrated by appropriate plots in the next subsections. 

3.1 Properties of the solutions corresponding to 

/ = y/2ao0^ 

3.1.1 General properties 

Figure 1 shows the typical behaviours of the solutions for 
discs with (dashed curves) and without (solid curves) a 
corona. Each curve is labeled by the appropriate values of 
(a ,l E /e). Behavior of the surface density is shown in the 
top left-hand plot of Figure 1. For the same input parame- 
ters, the presence of a corona enhances the surface density, 
and the surface density generally displays two regimes in 
the inner and outer regions. The size of the inner part is 
sensitive to the values of the input parameters. For a fixed 
viscosity coefficient, when the accretion rate increases the 
surface density in the inner part decreases, while the sur- 
face density at the outer parts increases in both cases with 
and without a corona. For a fixed accretion rate when the 
viscosity coefficient increases, the surface density decreases. 
Generally, we see that the corona becomes more effective 
when the viscosity and the accretion rate have larger values. 



Finally, at large radii the surface density is in proportion to 
R-°- 6 . 

The top right-hand plot of Figure 1 shows the ratio of 
disc thickness to the radius. The corona causes the disc to 
become thinner in comparison with the same disc without 
a corona. In the absence of a corona, for a fixed rate of 
accretion, viscosity increases the disc thickness. However, in 
a disc with a corona, so long as the accretion rate is constant, 
the thickness of the disc decreases in the inner part and 
increases in the outer part due to increase of the viscosity 
coefficient. With a fixed viscosity coefficient and increased 
accretion rate, the thickness of the disc increases in the inner 
part, but the thickness of the outer part decreases in both 
the case of a disc with a corona and the case without a 
corona. The disc thickness scales like R~ 13 at large radii. 

The bottom left-hand plot of Figure 1 shows that the 
presence of a corona serves to cool the disc because of the 
energy transported from the disc to the corona. Generally, 
for a fixed accretion rate when the viscosity coefficient in- 
creases, the disc temperature decreases because more energy 
is transported from the disc to the corona. But for a fixed 
viscosity coefficient when the accretion rate increases, then 
the temperature of the disc increases too. At large radii the 
temperature is changing in proportion to R~ 09 . 

The bottom right-hand plot of Figure 1 shows that the 
corona causes the disc to become more self-gravitating in 
the vertical direction in comparison with the case without 
a corona while, for a fixed accretion rate, when the viscos- 
ity coefficient increases, the self-gravity of the disc in the 
vertical direction decreases. However, for a fixed viscosity 
coefficient, when the accretion rate increases, there are two 
regimes, in both the case with and without a corona. At 
smaller radii, self-gravity of the disc in the vertical direction 
decreases, but as we go to the larger radii the vertical compo- 
nent of self-gravity becomes more significant. At larger radii 
the ratio £ varies in proportion to R 2 ' 75 . This plot shows 
that the vertical component of the self-gravity can be com- 
parable to the vertical component of the gravitational force 
of the central object at radii that are not large. Clearly, this 
radius is within a distance from the central object where the 
corona can exist. For instance for a pair (0.03, 1), the self- 
gravity of the disc is of the order of the gravitational force 
of the central mass at around R ~ 300i? s . For the (0.03, 10) 
case this happens at around R ~ 400R S . It clearly shows 
that self-gravity in the vertical direction is important for a 
disc-corona system even at not large radii from the central 
object where a corona exists. 



3.1.2 The effect of the self- gravity of the disc 

We study the possible effects of the self-gravity of the disc 
on the profiles of the physical variables of the disc in this 
subsection. Gravitational stability of the disc and the rule 
of the self-gravity is studied in Section 4 separately. Figure 
2 compares profiles of the solutions with self-gravity with 
those without self-gravity (KS solutions) . The input param- 
eters are similar to Figure 1 and each curve is labeled by 
a pair (ao,te/e). We see that the surface density at radii 
smaller than 10 3 Schwarzschild radius decreases due to the 
effect of self-gravity of the disc, irrespective of the viscosity 



© 0000 RAS, MNRAS 000, 000-000 



5 




Figure 1. Profiles of the physical variables for a disc without a corona (solid lines) and with a corona (dashed lines) vs. radial location 
in the disc with fi = 1.0, /i m = 0.6, k = 1 and Ms = 1. Each curve is labeled by the viscosity coefficient ao and the accretion rate l^/e as 
a pair (ooi'eA)- A different color is used for the curve corresponding to each pair of the input parameters. Here, the MN prescription 
for the dissipated energy into the corona has been considered. Since the amount of energy transported into the corona increases with 
the viscosity coefficient, the possible effects of the corona on the profiles are more evident for the higher viscosity coefficient. While the 
temperature of the disc decreases because of the corona, the surface density of the disc increases. 



coefficient and the accretion rate. Although the fraction of 
£ is not very large at these radii, the profile of the surface 
density changes significantly because of the self-gravity of 
the disc. 

Also, the self- gravity of the disc implies a thinner disc 
with a smaller opening angle. This reduction in the thickness 
of the disc is at all radii and for all the input parameters. 
In our analysis, we do not consider the vertical structure of 
the accretion disc. But in studies of the vertical structure 



of the accretion discs and analysis of the energy transport 
in this direction (e.g., via convection), one should note that 
the self-gravity of the disc significantly reduces the thickness 
of the disc according to our analytical solutions for a disc 
and corona system. However, the temperature of the disc is 
independent of the self-gravity of the disc for all the input 
parameters. However, the ratio of the gas pressure to the to- 
tal pressure, in particular for high accretion rates, increases 
because of the self-gravity of the disc. 
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Figure 2. Effects of the self-gravity of the disc on the profiles of the physical variables are shown in these plots by comparing our 
solutions with a corona (solid lines) and those solutions with a corona but neglecting the vertical component of the self-gravity of the 
disc (KS solutions, dashed lines). The input parameters are the same as Figure 1. Each curve is labeled by a pair (aa,l^/e). Different 
colors are used to denote the curves corresponding to the different pair of input parameters. Here, the MN prescription for the dissipated 
energy into the corona has been considered. We see reduction in the surface density of the disc at radii approximately smaller than 10 [i R s 
due to the self-gravity of the disc. Beyond this radius the surface density increases because of the self-gravity of the disc. However, the 
temperature of the disc is independent of the self-gravity of the disc. 



3.2 Properties of the solutions when / is a free 
parameter 

Svensson & Zdziarski (1994) presented analytical solutions 
for the disc-corona systems, in which a constant fraction / 
is dissipated into the corona. Now, we analysis properties of 
the physical variables of those solutions that correspond to 
a constant fraction /. Figure 3 shows profiles for /i = 1.0, 
jU m = 0.6, k = 1, oto = 0.3, Ze/c = 1 and Ms = 1. Each curve 



is labeled by the fraction /, which is assumed to be con- 
stant. For easier comparison, solutions without self-gravity 
are shown by dashed lines. Interestingly, properties of the 
solutions are qualitatively similar to the solutions with the 
MN prescription. Here, the surface density of the disc de- 
creases due to the self-gravity of the disc at radii smaller 
than 10 3 Schwarzschild radius. This reduction is indepen- 
dent of the amount of the energy transported to the corona. 
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Figure 3. Physical profiles of the disc when the fraction of dissipated energy into the corona is constant. The input parameters are 
H = 1.0, Hm = 0.6, k = 1, c*o = 0.3, ZeA = 10 and Mg = 1. Each curve is labeled by the fraction / (which is assumed to be constant) and 
is shown by different color. Solid lines are for the solutions with the vertical component of the self-gravity of the disc and the dashed lines 
show solutions when self-gravity of the disc is neglected. These curves are showing similar behaviors to when the fraction / is described 
by the MN prescription (Figures 1 and 2). Existence of the corona lead to an increase of the surface density and a cooler disc. Also, the 
surface density in the inner parts of the disc reduces because of the vertical component of the self-gravity of the disc, though the surface 
density increases beyond these regions. 



However, as the fraction / increases, the profile of the sur- 
face density shifts upward either with or without self-gravity 
of the disc. This means that as more energy is transported 
via magnetic tubes from the disc to the corona, the surface 
density increases. 

The thickness of the disc also depends on the self- 
gravity of the disc and the fraction /. The disc becomes 
thinner when more energy is transferred from the disc to the 



corona. Also, for a given fraction /, self-gravity of the disc 
causes the thickness of the disc to decrease. The thickness 
of the disc corresponding to the solutions with self-gravity 
is more sensitive to the fraction / compared with the solu- 
tions without self-gravity. Again, in this case, the existence 
of the corona implies a cooler disc, but the temperature of 
the disc is independent of the self-gravity of the disc. The 
ratio /3 also decreases at the inner parts of the disc due to 
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Figure 4. Profiles of the cooling time-scale t coo i (solid lines) and 
the Toomre parameter Q (dashed lines) vs. radial location in the 
disc. The curves are corresponding to the MN prescription and 
the input parameters are fi = 1.0, fi m = 0.6, k = 1, 7 = 2 and 
Mg = 1. Each curve is labeled by a pair (ao, Ze A) ar >d is showing 
by a color. 
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the self-gravity of the disc. However, this reduction is more 
evident for a smaller fraction /. 



4 GRAVITATIONAL STABILITY OF THE 
DISC-CORONA SYSTEM 

Having our analytical solutions for the steady-state struc- 
ture of an accretion disc with a corona, we can study the 
gravitational stability of such a system. Although most au- 
thors have been using the Toomre parameter as the main 
criterion in order to address the gravitational stability of ac- 
cretion discs (e.g., Goodman 2003; Goodman & Tan 2004; 
Nayakshin & Cuadra 2005), recent numerical simulations 
and physical considerations show that the cooling time-scale 
is another important physical factor (e.g., Gammie 2001; 
Lodato & Rice 2004). Also, the analysis of KS of the grav- 
itational stability of an accretion disc with corona is based 
on just the Toomre parameter (they had also neglected self- 
gravity of the disc). In this subsection, we study gravita- 
tional stability of our solutions considering the Toomre pa- 



Figure 5. Profiles of the cooling time-scale i coo i and the Toomre 
parameter Q vs. radial location in the disc for 00 = 0.03 (top) 
and Qo = 0.3 (bottom). Solid curves are corresponding to our 
solutions, but dashed lines are for solutions neglecting self-gravity 
of the disc (KS solutions). The curves are corresponding to the 
MN prescription and the input parameters are fj, = 1.0, /i m = 0.6, 
k = 1, 7 = 2 and Ms = 1. Curves corresponding to £eA = 10 are 
marked and the rest of the curves are for ZeA = !■ 



rameter and the cooling time-scale. Since solutions with free 
/ and the MN prescription show qualitatively similar behav- 
iors, the stability analysis is restricted to those correspond- 
ing to the solutions with the MN prescription. First, we cal- 
culate these parameters, then we do a parameter study. 



4.1 Toomre Parameter 

It was shown by Toomre (1964) that a rotating disc is sub- 
ject to gravitational instabilities when the Q parameter 
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Q 



si 



ttGS' 



(25) 



becomes smaller than a critical value, which is close to unity, 
where c s is the sound speed inside the accretion disc and Q = 
Q.K is the angular velocity. Having our analytical solution 
for the system, Toomre parameter of the model becomes 



Q = 4.945 x a 



7/10-3/10 



Ad, 



-13/10/fe . 

e ; 



-2/5 j-2/5^,-27/20 



x/3 (7^-24)/20 (1 _ v/2qo/3m/2) 3/10_ 



(26) 



This equation along with two algebraic equations (23) and 
(24) give Q as a function of the radial distance. Generally, 
Q is much greater than unity in the inner parts of the disc 
which implies these regions are gravitationally stable and do 
not fragment. But the Toomre parameter decreases with in- 
creasing radial distance so that Q reaches the critical value 
of unity at a self-gravitating radius which we denote it by 
R sg - Thus, all regions with 7? > R sg are gravitationally un- 
stable and may fragment. However, this is only true as long 
as the cooling time-scale also satisfies a criterion which will 
be explored in the next subsection. 

It can be shown that when the Toomre parameter is 
of order unity the ratio £ is around 4.8 in our model. With 
Q = 1, we have 



4-ivGp 



2c s Qk 
H ■ 



When combined with equation (2) this gives 



From (3) and c s = wp/ p, we have 



Cs 

Hfl k 



(27) 



(28) 



(29) 



Using equations (28) and (29), we obtain an algebraic 
quadratic equation for The only physical root of this 
equation is £ = 2(1 + \/2) ~ 4.8. This means that at a 
radius where the Toomre parameters is around unity, we 
have ( ~ 4.8, i.e., the force of the self-gravity of the disc in 
the vertical direction is 4.8 times the vertical component of 
the gravitational force of the central object. 

The most unstable wavelength for the Q ~ 1 disc is 
of the order of the disc vertical scale height H (Toomre 
1964). Thus, the most unstable mode has radial wave num- 
ber fe mu = (QH)^ 1 and so the mass of a fragment at R — R Bg 
becomes 



which gives 



"-mil 



(30) 



M frag = 0.126 x 10 34 a ( 1 / 5 A- 6 / 5 M- 4/5 (^)- 7 / 5 J- 7 / 5 r 3 - 6/5 



x/3 ( M -32)/io (1 _ p) 2 {1 _ ^ 2ao p^)- 6/5 . 



(31) 



4.2 Cooling time-scale 

There are actually different approaches to calculate the cool- 
ing time-scale. We can estimate the cooling time-scale r coo i 



using simple cooling laws as the ratio of the internal en- 
ergy density e = ECg/7(7 — 1) to the cooling rate A, i.e. 
Tcooi = e/A. Using equation (7), we can write 



A = (j^jhT^ = A^MJ(iJ). 



(32) 



On the other hand, angular momentum equation (6) gives 

a H{p^y /2 p (2 -^l 2 = ^-n K MJ(R). (33) 

Having equation (33), we can rewrite equation (32) as 

A = aomMp g a S ) M/ y 2 ^ )/2 . (34) 
Thus, our cooling time-scale becomes 

(35) 



Tcool = 1 1 /Pgas x(2- M )/2 

ao7(7 - 1) K p 



and non-dimensional cooling time-scale becomes 



tcool 



Tcool 



1 



-J 



(2-/0/2 



(36) 



to 3ao7(7 - 1) 

where 7 is the ratio of the specific heats and to = 3/fiK. 
According to Gammie (2001) and Lodato & Rice (2004) the 
fragmentation occurs when i coo i < 1. 

4.3 Gravitational instability: A parameter study 

When the disc cools very rapidly, the self-gravity of the 
disc does not have enough time to prevent the formation 
of bound objects in the disc. Since a corona extracts energy 
from the disc, we think, existence of the corona decreases the 
cooling time scale and the disc becomes more prone to frag- 
mentation. As we have already mentioned, in the KS anal- 
ysis of the gravitational instability of the disc and corona 
system, only the Toomre criteria has been used. Not only 
we consider this criterion of fragmentation, but the cooling 
time-scale condition is also considered in this paper. 

Figure 4 shows the cooling time scale and the Toomre 
parameter versus radial location in the disc. For the energy 
exchange from the disc to the corona, the MN prescription is 
used and the other input parameters are /1 = 1.0, /i m = 0.6, 
k = 1, 7 = 2 and Ms = 1. Approximately, when 7 = 5/3 
the cooling time-scale increases by a factor of two in com- 
parison to the case with 7 = 2. Each curve is labeled by 
a pair {aa,l^/e). Solid lines are for i coo i and dashed lines 
correspond to the parameter Q. In the inner parts of the 
disc, the Toomre parameter is very large, irrespective of the 
input parameters. Although the cooling time scale is less 
than unity (except for (0.03, 1)) in these regions, the Toomre 
condition for fragmentation does not satisfy and so the disc 
does not fragment. However, as we move toward larger radii, 
the Toomre parameter decreases down to the values smaller 
than unity. But the nondimensional cooling time-scale in- 
creases from small values at the inner parts of the disc to 
values greater or smaller than unity depending on the input 
parameters. Since in the MN prescription the amount of the 
transported energy from the disc to the corona is directly 
proportional to the viscosity coefficient ao, as this parame- 
ter increases, more energy is transported from the disc into 
the corona and the cooling time-scale decreases. According 
to Figure 4 when the viscosity coefficient is 0.03, the nondi- 
mensional cooling time-scale increases to values larger than 
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Figure 6. The self-gravitating radius R S g (in Schwarzschild radius R$) and the mass of the first clumps Mf rag (in solar mass) for fi = 1.0, 
/i m = 0.6, ft = 1, and Mg = 1. Top plots are corresponding to the solutions with MN prescription and bottom plots are for the solutions 
with fraction / as a free parameter. Solutions with self-gravity are showing by solid curves, but dashed lines are corresponding to the 
solutions without self-gravity. For solutions with free /, each curve is labeled by this parameter and we assume l^/e = 1. 



unity at regions of the disc where the Toomre parameter 
is less than unity, irrespective of the accretion rate. When 
accretion rate decreases, the increase of t coo \ becomes more 
smooth. However, for a larger viscosity coefficient 0.3 the 
cooling time-scale is lower than unity even at larger radii. 
So, in this case, the disc can fragment at regions where the 
Toomre parameter reaches to values lower than unity. 

Figure 5 shows possible effects of the self-gravity of the 
disc on the Toomre parameter and the cooling time-scale. 
We compare curves of Q and t C ooi for «o = 0.03 (top plot) 
and evo = 0.3 (bottom plot). While solid curves correspond 



to our solutions, the dashed lines are for solutions neglecting 
the self-gravity of the disc (KS solutions). As we discussed, 
our solutions for «o = 0.03 are gravitationally stable con- 
sidering cooling time-scale and Toomre conditions. However, 
we see that while the self-gravity of the disc slightly causes 
icooi to decrease in the inner parts, this time-scale increases 
at larger radii. Even though Q is less than unity, since the 
dimensionless cooling time is long, the disc will not collapse 
complete but will be in a self-regulated self-gravitational 
state. Also, the effect of the self-gravity is more evident at 
high accretion rates. But for qq = 0.3 the cooling time-scale 
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is smaller than unity at all radii. Again, the self-gravity of 
the disc slightly causes the cooling time-scale to decrease in 
the inner parts, and a longer time-scale in the outer parts. 
So, when the Toomre parameter becomes less than unity, 
the disc becomes gravitationally unstable. In this case, the 
self-gravitating radius decreases because of the self-gravity 
of the disc, in particular at the high accretion rates. 

Using equations (26) and (31) we can calculate the self- 
gravitating radius R Bg at which the Toomre parameter be- 
comes unity and the mass of the fragments Mf rag at this 
radius. Figure 6 shows results of the calculations for a sys- 
tem with a central mass of 1O 8 M0, both for the solutions 
with MN prescription (top plots) and solutions with free 
fraction / (bottom plots). The other input parameters for 
the MN prescription are /i = 1.0, /i m = 0.6, k = 1 and 
Ze/e = 1 and 10. For comparison, while solid curves corre- 
spond to our solutions, the solutions neglecting self-gravity 
of the disc are shown by dashed curves. For the viscosity 
coefficient ao, we restrict to those values which correspond 
to a nondimensional cooling time-scale lower than unity. For 
values of the coefficient ao close to 0.5, the self-gravitating 
radius decreases. But the location of this radius is more sen- 
sitive to the viscosity coefficient for low accretion rates either 
with self- gravity or without it. Interestingly, the self-gravity 
of the disc causes the self-gravitating radius to decrease, in 
particular for low accretion rates. But as the viscosity coeffi- 
cient increases, the self-gravitating radii of the solutions with 
self-gravity and without self-gravity become closer. On the 
other hand, we note that in the MN prescription the fraction 
of dissipated energy into corona is directly proportional to 
Qo- Thus, as more energy is transported from the disc to the 
corona, the location of the self-gravitating radius becomes 
more independent of the self-gravity of the disc. Figure 6 
(top; right plot) shows the mass of the fragments at the self- 
gravitating radius in the case of the MN prescription. The 
mass of the fragments increases significantly when the ac- 
cretion rate is high and, generally, the mass of the fragments 
increases because of the existence of a corona (KS). But the 
self-gravitating radius decreases due to the existence of the 
corona either with low or high accretion rate. Consequently, 
the presence of a corona makes the disc more gravitationally 
unstable (KS). Also, the self-gravity of the disc significantly 
lead in reduction of the mass of the fragments for all values 
of the viscosity coefficient or accretion rate. 

For comparison, the bottom plots of Figure 6 show the 
self-gravitating radius and the mass of the fragments for 
the solutions with free parameter /. Each curve is labeled 
by a corresponding fraction /. The other input parameters 
are n — 1.0, pi m = 0.6, k = 1 and Ze/e = 1. Again, so- 
lutions without self-gravity of the disc are shown by the 
dashed curves. For the effect of the self-gravity, in this case, 
behavior of the variations of R sg and Mf rag with the input 
parameter are qualitatively similar to the solutions with MN 
prescription: reduction of the self-gravitating radius due to 
the self-gravity of the disc, but as more energy is transported 
from the disc to the corona, location of the self-gravitating 
radius becomes independent of the self-gravity of the disc. 
However, the mass of the fragments significantly decreases 
due to the self-gravity of the disc, irrespective of the input 
parameters, notably the parameter /. 



5 DISCUSSION 

Generally and for simplicity, the radial or the vertical com- 
ponents of the gravitational force due to the self-gravity of 
the disc have been neglected in most of the analytical stud- 
ies of self-gravitating accretion discs. However, some authors 
emphasize on the importance of the gravitational force due 
to the self-gravity of the disc (e.g., Paczynski 1978; Hure 
1998; Bertin & Lodato 1999). To our knowledge, most of the 
previous studies of star formation in AGNs or the Galactic 
center do not consider the components of the gravitational 
force due to the self-gravity of the disc (e.g., Goodman & 
Tan 2004; Nayakshin & Cuadra 2005). Although KS studied 
the possible effects of a corona on the gravitational stabil- 
ity of an accretion disc, they also neglected the components 
of the self-gravity of the disc itself. In this work, we pre- 
sented sets of analytical solutions for the steady-state struc- 
ture of self-gravitating accretion discs with corona. We also 
considered the vertical component of the gravitational force 
due to the self-gravity of the disc. The energy exchange be- 
tween the disc and the corona has been considered either as 
a function of the physical variables of the disc, or as a free 
parameter of the model. Although there are possible mass 
and angular momentum exchanges between the disc and the 
corona, we have neglected these effects just for simplicity in 
order to present a self-consistent model for a disc-corona sys- 
tem. Having the analytical solutions, we studied properties 
of the solutions, more importantly gravitational stability of 
the disc. For the conditions of fragmentation, not only the 
Toomre parameter but the cooling time-scale are considered 
in our analysis. 

We showed that because of the energy exchange, the 
corona lead in a cooler the disc for all the input parameters. 
This reduction of the temperature due to existence of the 
corona, decreases the thickness of the disc. Since there is 
no mass exchange between the disc and corona, the surface 
density increases and as a result, the vertical component of 
the gravitational force due to the self-gravity of the disc be- 
comes more effective. When the Toomre parameter reaches 
to around unity, the ratio £ is around 4.8 which means that 
the vertical component of the gravitational force due to the 
self-gravity of the disc is a few times larger than the vertical 
component of the gravitational force of the central object 
at the self-gravitating radius. Thus, this term may modify 
properties of the solutions as our solutions show. For exam- 
ple, self-gravity of the disc in a disc-corona system causes 
the surface density to decrease at the inner parts of the disc, 
though the temperature of the disc is unchanged when the 
self- gravity of the disc is considered. 

Considering Toomre parameter and cooling time-scale 
conditions for fragmentation, solutions with small viscos- 
ity coefficients are stable to fragmentation because in the 
inner parts of the disc Q > 1 and when Toomre parame- 
ter reaches to unity, the cooling time-scale becomes greater 
than one. This result is independent of the accretion rate 
or other input parameters. But for large viscosity coefficient 
the cooling-time scale is smaller than one for all radii, and 
so, once the Toomre parameter becomes less than one the 
disc may fragment. However, the cooling time-scale slightly 
decreases at the inner parts of the disc because of the self- 
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gravity of the disc and for high accretion rates, this time- 
scale increases. 

Also, self-gravitating radius and the mass of the frag- 
ments at self-gravitating radius decrease because of the self- 
gravity of the disc. However, change of the self-gravitating 
radius due to the self-gravity is negligible when the amount 
of the energy which is transported from the disc to the 
corona is increased. But reduction of the mass of the frag- 
ments because of the self-gravity is independent of the frac- 
tion /. These results are obtained either with free / or MN 
prescription for the fraction of energy exchange. Although 
we estimate the mass of the fragments in an accretion disc 
with corona around supermassive blackholes, we can not de- 
termine whether these fragments are sufficiently long-lived 
to be able to evolve in protostars over many orbital periods. 
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